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The out-of-plane current-voltage (I-V) characteristics of Bi2212 are studied in experimental envi- 
ronments of different heat transfer efficiency, allowing practical separation of intrinsic and extrinsic 
phenomena. Intrinsic (heating-free) response is Ohmic in the normal state of Bi2212, while its 
resistance, R=V/I, is found to be a good practical measure of the mean temperature of the sample 
in the overheated case. A self-heating model proposed for the latter case provides a qualitative and 
quantitative description of key findings of intrinsic tunnelling spectroscopy including (pseudo)gaps, 
quasiparticle and normal state resistances. The model also naturally explains the 'superconducting' 
gap closure well below T c of the material as well as its survival at a magnetic field significantly ex- 
ceeding H C 2- The generic shape of the individual branches of the brush-like part of I-V established 
under conditions of negligible overheating suggests a phase-slip origin of the intrinsic Josephson 
effect (IJE). 

PACS numbers: 74.45.+C, 74.50.+r, 74.72.-h, 74.25.Fy, 



I. INTRODUCTION 

Despite significant efforts to understand the mecha- 
nism of high temperature superconductivity (HTSC), no 
consensus has been achieved either in theory [l| or ex- 
perimentally. In particular, conflicting conclusions on 
the value and temperature dependences of the supercon- 
ducting gap, A s , and the normal state pseugo-gap, A p , 
in the electronic density of states (DOS) have been drawn 
from conventional spectroscopy [E El IE IE ancl 
from 'intrinsic' tunnelling (c/ ; 8] and references therein). 
Intrinsic tunnelling (ITs) is believed to occur in highly 
anisotropic superconductors such as Bi2212 and Bi2201, 
manifesting itself in the peculiar shape of an out-of-plane 
I-V (IVC). The low- bias brush-like part of such I-V is 
highly hysteresial and becomes single-valued at higher 
voltages, where gap-like peculiarities are located. In re- 
semblance of 'conventional' Josephson junctions' I-V, in- 
dividual branches of that brush were found to be highly 
nonlinear. 

This similarity has led to a belief that layered HTSC 
such as Bi2212 represent natural stacks of atomic scale 
intrinsic superconductor- insulator-superconductor (SIS) 
- Josephson junctions. Correspondingly, it is generally 
believed that the c-axis (out-of-plane) charge transport 
in layered HTSC occurs via a sequential tunnelling of 
electrons or Cooper pairs between Cu-0 planes. Some 
important conclusions have been drawn from results ob- 
tained with this method, where Josephson coupling be- 
tween Cu-0 planes was assumed to be the only cause 
of the peculiar shape of I-V. This assumption, however, 
is not beyond dispute from the experimental viewpoint 
[E ^E El- I n particular, unlike in conventional tun- 
nelling, the characteristic features of ITs are seen at high 
heat loads, P, (estimated as VI/ A, where A is the area 
of the sample) in excess of kilowatts per cm 2 at T^4.5K. 



This value significantly surpasses the critical load for liq- 
uid 4 He, ~lW/cm 2 , thus signalling possible heating, es- 
pecially since the poor thermal conductivity of all HTSC 
makes them particularly prone to local overheating. Nev- 
ertheless, heating issues in IT studies of HTSC were 
mostly misinterpreted or ignored until recently, when di- 
rect experimental evidence of the heating origin of I-V 
non-linearities in Bi2212 mesas was obtained [l0| . 

The majority of IT investigations use bridge- like struc- 
tures (mesas) fabricated from HTSC crystal or film. Un- 
like tunnel junctions, these mesas, like any superconduct- 
ing weak links, are characterised by a high current den- 
sity, which generates a nonequilibrium state in the struc- 
ture [12J. The stationary and non-stationary properties 
of the nonequilibrium state in conventional superconduc- 
tors have been a subject of systematic investigation. In 
particular, phase-slip centers (PSC) of different dimen- 
sionality have been studied in great detail 0, Q, 
as have a rich variety of nonlinear phenomena caused by 
Joule self-heating. The latter often mimics Josephson- 
like behaviour |l(j| so that care must be taken to distin- 
guish intrinsic properties from extrinsic. 

Below I report on an experimental study of c-axis I- 
V of Bi2212 single crystals and mesa-structures of the 
lateral in-plane size, W=10-300/ito, performed both un- 
der conditions of negligible self-heating and in the oppo- 
site scenario, where I-V characteristics are dominated by 
Joule self-heating. A self- heating model of IT is proposed 
on these experimental grounds and is found to provide a 
quantitative description of IT phenomenon. Examples 
drawn from the most accurate investigations will show 
that the majority of results obtained with the IT tech- 
nique could be reproduced using the temperature depen- 
dence of sample resistance in the normal state, i?/v(T), 
and with the assumption that the heating is the only 
cause of non-linearities observed in current-voltage char- 
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acteristics. I will also briefly discuss heating-free data 
and possible causes of multiple branching. The technique 
employed for study of the quasi-stationary I-V, together 
with some other experimental details, can be found else- 
where [n|. 



II. RESULTS AND DISCUSSION 

Our discussion begins with the experimental out-of- 
plane I-V measured under conditions of both negligible 
self-heating and in the opposite scenario, where Joule 
self-heating dominates. For the latter case, the self- 
heating model is proposed. A few examples drawn 
from reliable publications will then demonstrate that this 
model self-consistently explains the most significant find- 
ings obtained with IT spectroscopy, including the pseu- 
dogap, the superconducting gap and their coexistence, as 
well as the contrasting effect of magnetic field on A s and 
A pg . The following sections consider low-bias IT results 
and some smoking guns for the model; and the last sec- 
tion is devoted to the 'heating-free' I-V and the possible 
cause of the brush-like type of out-of-plane I-V in Bi2212. 



A. Experimental basis of the model 
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FIG. 1: Inset A shows the typical I-V in the presence of heat- 
ing. The symbols in B present the same I-V, re-plotted as V/I 
vs W, together with ac-R(T) of the same sample measured at 
low bias (solid line, the top x-scale). C Heat loads required 
to heat the crystal from Tb to 100K. 



At first, we considered how the heating-free I-V would 
look. In our studies, we followed several simple guide- 
lines. Any I-V is inevitably accompanied by Joule dissi- 
pation which can escape only to the bath. Mean over- 
heating therefore depends on the Joule heat and the ef- 
ficiency of the heat transfer into the bath, which can be 
modified significantly. However good the thermal con- 
tact, special attention must be paid to the initial (low- 
bias) part of the characteristics, where overheating is 
smaller. 

Indeed, the measurements performed at sufficiently 
small heat loads in the normal state of the Bi2212 crystal 
were found to be virtually unaffected by self-heating. The 
linearity of the out-of-plane I-V observed in this case over 
4-9 orders of heat load magnitude and reasonable con- 
cordance between the threshold values P c , at which the 
I-V linearity breaks down and where the Joule heating 
of the crystal becomes measurable, suggest the intrinsic 
(heating-free) response to be Ohmic in the normal state 
of Bi2212. 

This approach is not directly applicable below T c since 
the superconducting critical current, I c , imposes a lower 
limit on the measuring current. A non- vanishing I c of- 
ten causes transition of the superconductor into a self- 
sustaining regime in which the Joule heat release becomes 
sufficient to raise the temperature above T c ^(|. For 
this reason, below T c the 'heating-free' I-V can only be 
addressed reliably under conditions of sufficient I c sup- 
pression. Both temperature and magnetic field make it 
possible to suppress / c , thereby providing a practical tool 
to distinguish intrinsic phenomena from self-heating. 



In agreement with the earlier findings [ToL IhH | , it was 
found that at sufficiently low temperatures self-heating, 
which may occur at currents much lower than I c , de- 
stroys the superconducting state in Bi2212. Moreover, 
my experiments demonstrate that I-V similarity with 
Josephson or tunnel characteristics does not necessarily 
imply them to be of the same origin. For example, Joule 
overheating of the crystal is largely responsible for the 
Josephson-like shape of I-V in Fig.lA. 

In terms of heating issues, it is appropriate to con- 
sider R=V/I as a function of the heat load, P = VI/ A, 
rather than I-V only (see for details). The shape 
of the I-V curves accordingly re-plotted recalls that of 
the experimental R(T) of the same sample. Moreover, 
these different dependencies could be superimposed by a 
simple re-scaling, as shown in Fig. IB, which presents the 
same I-V properly re-plotted together with ac-R(T) of 
the same sample. Reasonable quantitative concordance 
between the measured and calculated R(T) is evident 
from this figure. This finding suggests that, where the 
crystals' overheating is determined by the heat load, the 
sample resistance, R=V/I, is a good practical measure of 
the overheating. 

Moreover, this finding allows an experimental estimate 
of the characteristic heat loads required to heat the crys- 
tal from the arbitrary bath temperature, Tg, to a certain 
target temperature, Tt- For example, Fig.lC shows the 
experimental dependencies obtained for Tt = lOOA' on 
the same sample in different experimental environments. 
These curves obey a linear fit (broken lines in Fig.lC), 
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whose slopes depend on the heat transfer efficiency of the 
set up used. 

Fig.lC provides experimental evidence of the appli- 
cability of Newton's Law of Cooling (1701) to layered 
HTSC. This result, together with the generic out-of- 
planc I-V linearity in the normal state of the layered 
HTSC, are of principal importance for the self-heating 
model [Toj , which assumes that heating is the only cause 
of non-linearities observed in current- voltage characteris- 
tics. Below it will be shown that the majority of results 
obtained with IT technique could be reproduced using a 
simplified linear version of the model. 



B. The pseudogap, A pg 

Let us first address the origin of the pseudo-gap, A p , 
observed with IT spectroscopy (see, for example, [T^ and 
references therein). From a typical example it may be 
seen that the 'tunnelling characteristics' measured in |8( 
over a wide temperature range in the normal state of 
layered HTSC could be reproduced quantitatively using 
experimental normal state resistance R(T) of the same 
sample (reproduced in the insert to Fig. 2), with the as- 
sumption that the mesa heating caused by Joule dissipa- 
tion is the only reason for effects observed at high bias 
(see also H3). Using Newton's Law of Cooling (1701), 
the temperature of a thin mesa is given by 

T = T B + IV/(Ah), (1) 

where Tb is the temperature of the coolant medium (liq- 
uid or gas) and h is the heat transfer coefficient. Ac- 
cording to Eq.(l), monitoring I-V at a given bath tem- 
perature Tb results in a sample temperature rise that 
entails non-linearity in I—V/R(T), which is gap- like if 
dR/dT<0. The alL/dV curves thus constructed resemble 
those presented in fig. 4 of Q and reveal a quantitatively 
similar variation of conductance. This similarity allows 
for an estimate of the heat transfer (hAjc^WfxW / 'K and 
overheating at high bias, ^80K for Tb=20QK ; the set of 
curves accounting for this coefficient is shown in Fig. 2. As 
that Figure clearly shows, the 'heating' spectra taken at 
Tb<T* reveal a 'pseudo-gap' which disappears entirely as 
soon as Tb>T*. Interestingly, the IVCs in Fig. 2 are non- 
linear even above T*. In perfect agreement with the ex- 
periment, the heating spectra have an inverted parabolic 
shape, which was originally attributed to the presence of 
van-Hove singularity close to the Fermi level, cf. eg. [l9j . 

Thus, Eq.(l) provides a natural and adequate explana- 
tion of the puzzling findings of Q . Accordingly, the non- 
linear I-V characteristics observed by 8] in Bi2201 are 
undoubtedly related to the temperature dependence of 
the normal state c-axis resistance of their sample rather 
than to a (pseudo)gap in the tunnelling density of states 
(DOS). Moreover, there is little doubt that R(T) modifi- 
cation caused by doping is responsible for the rich zoology 
of the 'pseudogaps' observed in Chalmers p4| . 
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FIG. 2: dl/dV obtained from R N (T) of [|] (shown in the 
insert) assuming Joule heating origin of I-V nonlinearities. 



C. A s and A pg coexistence 

The validity of Eq.l down to 30K, verified in per- 
mits its application to the results originally attributed to 
the superconducting gap. A lack of information about 
normal state resistance in the absence of superconduc- 
tivity, i?j\r(T) makes this task more problematic. Nev- 
ertheless, as was shown elsewhere [l£|, the same model 
provides reasonable agreement with IT spectra mea- 
sured below T c using the experimentally verified assump- 
tion 0, that the 'sub-gap' IT-resistance, R sg = 
dV/dL\p^ , obtained in zero field is a good practical 
measure of the Rn of a single sample in conditions where 
its superconductivity is destroyed, for example by suf- 
ficiently high magnetic fields. The resulting 'heat ing ' 
dl/dV curves resemble the experimental ones from |20j 
and allow estimation of both the 'gap' and the heat trans- 
fer, (hA) « 25 — 33fiW/K. The 'heating' spectra reveal a 
gap-like peculiarity whose temperature dependence is re- 
markably similar to that in |2(j. Hence, Eq.(l) provides a 
natural explanation of some of the IT-results attributed 
to the superconducting state pcj . 

However, the 'heating' IVC reveal one set of gap-like 
or S-shaped peculiarities, provided that Rn(T) is a well 
behaved function as observed 0, l^. Although this 
agrees perfectly with the majority of IT results, there 
are a few works which claim that their IT-spectra 'evi- 
dence the coexistence of the superconducting gap and the 
pseudogap' [l^. |24|. At first glance, these observations 
cast doubts on the model's applicability. More close ex- 
amination, however, reveals that the gaps' coexistence is 
observed in 3-probe measurements where contact resis- 
tance, i? c , inevitably contributes to the net R(T). It is a 
priori evident that R C {T) may reveal vastly different be- 
haviour depending on the type of particular contact; see 
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for example [25|. Besides, R C (T) is inevitably affected 
by development of the true superconducting gap of the 
material. Moreover, for certain types of contacts the con- 
tact resistance may become highly non-Ohmic below T c 
26] . For these reasons, the lack of R C (I, T) curves 
increases the ambiguity in low temperature behaviour of 
the net R(T). 

Although the influence of the 3-point configuration on 
the 'heating' IT-spectra might be evaluated for a general 
case of NS contact of an arbitrary barrier strength using 
IVC by [2{| together with the normal state R(T) mea- 
sured in the absence of superconductivity 0, 153, l23| . 
this task falls beyond the scope of this article and will be 
addressed elsewhere. 

Below I will outline on a qualitative level how the R c 
modification caused by the superconducting gap (SG) de- 
velopment might influence the resulting IT spectra. I 
will show through a simple example that this R c modi- 
fication not only affects the net 3-point resistance R(T), 
but modifies the 'heating' spectra radically. Neglect- 
ing R c nonlinearity for the sake of simplicity, let us as- 
sume that the SG contribution to the total contact re- 
sistance, i? c , in the T c vicinity might be modelled as 
5R C (T) =SR C (0)(1-T/T C ). This contribution may cause 
little effect on the net R(T) provided that 5R C (0) is suf- 
ficiently small. To visualize this statement, 3-point R(T) 
curves (broken lines in Fig. 3) are constructed from the 
well-behaved 'mother' one (shown by the solid line) for 
a representative set of SR c (0). For the sake of trans- 
parency, this construction uses the same 'mother' R(T) 
as in Fig. 2 and assumes T c = 60K, as shown by the ar- 
row in Fig. 3. At first sight the thus distorted curves seem 
to differ insignificantly from the undisturbed one even 
for the highest SR c (0) = 30f2 among those presented in 
Fig. 3. Unlike R(T), the heating spectra are drastically 
affected by this modification. As is clearly seen in the in- 
sert to Fig. 3, unlike the 'mother' heating spectrum, the 
spectra modelled for Tb < T c reveal both a 'pseudo-gap', 
V pg , and a 'superconducting' gap, V s . 

Thus, this figure reproduces the most significant find- 
ing of [l3)l23 : it shows that a sharp 'superconducting' 
peak at V s emerges on top of the hump at V pg associ- 
ated with the normal state pseudogap in [l9|. The lat- 
ter feature is evidently a robust signature of the normal 
state R(T), while the 'superconducting' peak is prob- 
ably caused by the 3-point configuration employed in 
pi fillip . The extrinsic origin of the V s peak in Fig. 3 
is additionally illustrated by the strong dependence of 
its amplitude on 5R C (0) value. Moreover, more sophis- 
ticated analysis has shown the amplitude and position 
of this peak V s , to be very sensitive to the particular 
shape of SR C (T) as well as to R c nonlinearities j^. So 
our model provides a natural explanation of the puzzling 
findings of 0, 0, 0]. However, given the fundamental 
importance of the conclusions of 0], further evidence 
of the extrinsic origin of the A s and A pg coexistence is 
desirable. 

Fortunately, one study, [24|, allows an estimate of the 
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FIG. 3: Modelling the 3-probe configuration induced coexis- 
tence of A s and A P9 within the heating scenario. The solid 
line represents the 'mother' well-behaved R(T) from the inset 
to Fig.2. Broken lines are constructed from this R(T) assum- 
ing that a superconducting gap development below T c = 60K 
(shown by the arrow) contributes to the contact resistance as 
SR C (T) = SR C (0)(1 - T/T c ) with 5R C (0) = 6; 12; 18; 24; 30fi 
respectively. Insert shows dl/dV obtained from the corre- 
sponding R(T) curves in the main panel assuming Joule heat- 
ing origin of I-V nonlinearities. 



net resistance of the bridge, thus offering an independent 
test of our conclusions. Unlike in the case of direct mea- 
surements of normal state resistance in the absence of 
superconductivity, eg. |lOj . and in qualitative agreement 
with the a priori expectations mentioned above and il- 
lustrated in Fig. 3, the resulting 3-point R(T) from [24[ 
reveal a kink-like change below T c and a temperature- 
dependent discrepancy between the normal state R(T) 
extrapolation towards T < T c and the low temperature 
estimates by [24] . In resemblance of the calculated curves 
in Fig. 3, the experimental discrepancy increases rapidly 
as the temperature lowers further below T c . Remarkably, 
the relative magnitude of this discrepancy exceeds that in 
Fig. 3 by more than one order, thus additionally support- 
ing our assumptions. Moreover, the spectra constructed 
with Eq.(l) from the net R(T) by [24J are qualitatively 
similar to those presented in the insert to Fig. 3. In re- 
markable resemblance of the experimental IVC by the 
Chalmers group [24(, these spectra reveal both a 'pseudo- 
gap', Vpg, and a 'superconducting' gap, thu s ag ain repro- 
ducing the most significant finding of [l^ . l24j : a sharp 
'superconducting' peak at V s emerges on top of the h ump 
at Vpg associated with the normal state pseudogap in . 

Unlike the computed curves in Fig. 3, the low tempera- 
ture measurements by [24| reveal noticeable scatter, thus 
permitting slightly different R(T) approximations. Inter- 
estingly, in perfect agreement with |27j . vastly different 
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amplitudes of the 'superconducting' peak and its posi- 
tions, V s , were obtained with Eq.(l) from these R(T) 
approximations, which are compatible with the error bar 
of |2J|. Unlike the Chalmers group, I believe that the V pg 
feature is a robust signature of the normal state R{T), 
while the 'superconducting'peak is caused by the 3-point 
configuration employed in 0, E, 0] • 

Thus, Fig. 3 together with the 'heatin g' s pectra cal- 
culated from the net 3-point R(T) from [2J] provide a 
natural explanation of the puzzling findings of [24| . In- 
terestingly, the heat transfer efficiency (hA)~8AfiW/ K 
and overheating at high bias, >200K for Tb=60K esti- 
mated for the case of [24j correlate reasonably with the 
estimates obtained in the previous section. Moreover, 
complementary to the 'heating' spectra, our model al- 
lows reconstruction of R(T) dependence from intrinsic 
tunnelling I-V obtained at fixed Tb- Remarkably, the 
same hA allows perfect superposition of the R(T) thereby 
constructed with Eq.(l) from I-V by |24j with experimen- 
tal R(T) of the same sample. Thus, the results of this 
section demonstrate the self-consistency of the model. 

Hence, it is most probably fallacious to identify the 
symmetric peaks in IT-spectra with the superconducting 
energy gap A s or normal state pseudo-gap A p in the elec- 
tronic density of states (DOS). Indeed, experimental IT 
spectra undoubtedly originate from Joule heating and are 
related to the peculiar temperature dependence of the net 
resistance of the bridge. The three-point configuration 
developed in Chalmers is most likely responsible for the 
resulting peculiar if not pathological R(T). The 'heating' 
spectra might of course also be affected by the contacts' 
I-V deviations from linearity. However, this contribution 
cannot be verified reliably because of the lack of contact 
characterisation in |24| . 



D. Magnetic field influence on A s and A pg 

Let us turn now to the puzzling effect of the magnetic 
field on the IT features ascribed to the superconduct- 
ing gap, V^, and the c-axis pseudogap, V pg . Recently, 
Krasnov [l9( observed that V s and V pg have different 
magnetic field dependencies. Ignoring self-heating, Kras- 
nov considered these observations to represent a crucial 
test for the superconducting origin of the two gaps in 
DOS of HTSC. However, the thermal load estimate even 
for HgBr2-Bi2212 crystal at Helium temperature exceeds 
KILOWATTS/cm 2 , so that this assumption is probably 
incorrect. 

Moreover, I will show below that the same model al- 
lows a natural explanation of this IT finding, namely the 
contrasting effect of the magnetic field, B, on V s and V pg , 
where the first characteristic was found to be much more 
sensitive to the field 0,12^. In my opinion, the normal 
state magnetoresistance (MR) of the substance investi- 
gated and its temperature dependence |22j are responsi- 
ble for that finding. Indeed, as shown above, to reach 
V s and V pg features, one has to warm the sample almost 
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FIG. 4: dl/dV obtained from R(B, T) of @ (shown in the 
insert) assuming Joule heating origin of I-V nonlinearities. 

up to T c and T*, the minimum in R(T), respectively. 
The negative normal state MR in Bi22I2 is experimen- 
tally found to be small at T ~ T* [23, HI] , thus ensuring 
negligible V pg (B ) d ependence, in perfect agreement with 
the experiment |l9L l28j ] . However, the magnitude of the 
normal state MR grows rapidly on temperature lowering 
[23,|2^|, so that the noticeable field influence on V s i s also 
explained naturally within the heating scenario [lfllril |. 

Fortunately, the data of [23 allow a cross-check of 
the model predictions, since 3-point R(T) curves are 
presented for a set of fields. Of course, the unknown 
R C (B, T, I) in 28] prevents any reliable MR estimate for 
the HTSC investigated. However, these curves allow an 
estimate of the net 3-point MR which is responsible for 
the magnetic field influence on the IVC. This MR is huge, 
30%(@14Tesla) at 4K and it drops below the resolution 
of |28j as the temperature increases towards T c , as seen in 
the insert to Fig. 4. Evidently, these R(T) curves allow es- 
timation of the effect of magnetic field on the 'heating' IT 
spectra. Fig. 4 shows dl/dV spectra obtained with Eq.(I) 
from R(T) measured at two fields. Qualitative similarity 
with the results of 0,H^] is evident from Fig. 4. Indeed, 
the heating spectra in Fig. 4 evidence gap suppression by 
magnetic field and temperature in reasonable agreement 
with the experiment [l^ . |28[ . I believe this similarity 
favours the self-heating origin of the results obtained in 
Chalmers especially given the unknown field, current and 
temperature dependence of their contact resistance, R c . 

Interestingly, the heat transfer efficiency 
(hA)~8.5fj,W/ K estimated for this experiment is 
strikingly similar to other cases considered here, which 
provides extra support for the model. So the most 
significant findings of high-bias ITs could be broadly 
understood within the heating scenario. 
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In conclusion, temperature dependent normal state 
magnetoresistance is the most likely cause of the con- 
trasting response to magnetic field of the two gaps ob- 
served in intrinsic tunneling experiments |l9l |28[ . This 
provides a natural explanation of the V s closure well be- 
low T c of the material as well as its survival at magnetic 
fields which significantly exceed H r 2, the upper critical 
field of the material [I |H IH ^ ^3 Although the 
low-bias results reported by this group might be less af- 
fected by self-heating 0, On], these are of little rele- 
vance because of the 3-point configuration and unknown 
field-temperature-current dependence of the contact re- 
sistance employed in 0, 0, |28| . 



E. Low-bias IT data 

Undoubtedly, there is a certain level of dissipation at 
which intrinsic phenomena may dominate heating. Ex- 
tra experiments are required in order to address this is- 
sue unequivocally. However, certain preliminary conclu- 
sions may be drawn from analysis of the so-called 'low- 
bias' IT-curves. To simplify the task, let us consider 
IT results obtained on the Dayem-bridge-like HTSC- 
structures. These structures benefit from spatial sepa- 
ration between the junction and the contacts, since the 
Joule heating of the contacts might influence the mean 
temperature of the junction less than in the case of the 
stacked geometry of mesas considered in the previous sec- 
tions. Thus, these objects are highly suitable for the pur- 
pose of this section. 

Similar to other 4-point IT-studies, the current-voltage 
characteristics of these bridges reveal one gap-like pe- 
culiarity, V s and exceptionally well defined and nearly 
temperature-independent Ohmic resistance at very high 
bias, which was originally attributed to the normal- 
state tunnel resistance, R n . However, as was first noted 
in 0, the very high levels of dissipation required to 
achieve these characteristic features at helium tempera- 
tures, namely, 7.3 and 210kW/cm 2 respectively, evidence 
the heating origin of these features. Fortunately, com- 
plementary to the high-bias IVC, the authors of [23 a l so 
measured a low-bias I-V (see Fig.5A) which may be less 
affected by heating. It will be shown below that our 
model allows discrimination between regimes of strong 
self-heating and those where intrinsic phenomena may 
dominate heating. It will be seen that our model al- 
lows estimation of the highest (critical) heat load that 
this particular sample can bear at Tb without notice- 
able overheating. Thus, it allows discrimination between 
regimes of strong self-heating and those where intrinsic 
phenomena may dominate heating. 

Evidently, monitoring I-V at a given base tempera- 
ture, Tb, results in a sample temperature rise as soon as 
P exceeds the critical heat load P c . Indeed, as is seen 
from Fig.5B, there is a well defined threshold level, P c , 
of the heat load below which R(P) is flat, while it drops 
rapidly at P > P c . The R(P) curves in the latter case 
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FIG. 5: A: low-bias I-V of the bridge #2 from Fig. lb and 3 
of |3^ | are shown by the solid and broken lines respectively; 
B: heat load dependence of bridge #2 resistance, R=V/I, 
obtained from Fig.5A. 



are probably caused by Joule self-heating. Indeed, simi- 
larly to other IT examples considered here, these curves 
resemble R(T) by |32] and reveal a quantitatively similar 
variation of resistance, which allows an estimate of the 
heat transfer coefficient /i^30-60Wcm -2 K -1 (the scat- 
ter reflects the discrepancy between Fig. lb and 3 in |3^| 
which is also clearly seen in Fig. 5 A here). The set of 
curves accounting for this coefficient correlate reason- 
ably with the measured R(T) curves, thus confirming 
the heating origin of the falling part of R(P) dependence 
in Fig.5B. Thus, Eq.(l) provides a natural explanation 
even of the 'low-bias' findings of [32 and suggests that 
the 'heating-free' I-V is linear, in agreement with the con- 
clusions drawn from the analysis of IVC measured over 
an extended range [l8j . 

To conclude, the character of the Joule heating of a 
sample with typical c-axis R(T) dependence is respon- 
sible not only for V s and R n but also for the majority 
of the low-bias data by |3^. Moreover, our model allow 
evaluation of the critical heat load below which intrin- 
sic phenomena may dominate heating. The linearity of 
the thereby determined intrinsic (heating-free) response 
of Bi2212 from [33 suggested by Fig.5B has little in com- 
mon with the Josephson behaviour taken for granted by 
IT advocates. 



F. Some more smoking guns for the model 

The self-heating model shows that the vast variety 
of results obtained with intrinsic tunnelling spectroscopy 
most certainly originate from Joule heating of the sample 
and are related to the temperature dependence of normal 
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state resistance. The remarkably similar heat transfer co- 
efficients estimated here for various Bi2212 and Bi2201 
microbridges explored by different groups worldwide pro- 
vide additional support for the model. Moreover, there 
are several more 'smoking guns' which demonstrate that 
the model is correct. 



In particular, Eq.(l) suggests that, providing h is con- 
stant, mesas with different areas would require the same 
heat loads to achieve the same I-V features. Indeed, strik- 
ingly similar loads are required to achieve the character- 
istic features of the I-V of Bi22I2 mesas of a vastly differ- 
ent area, 25-80/im 2 , at the same base temperature, 4.3K 
namely, 3.3kW/cm 2 and 18. 5kW/ cm 2 for those those as- 
cribed in [24| to the superconducting gap, V s , and to the 
normal state Ohmic response, i?„, respectively. 

Magnetic field offers one more test of the model, since 
it modifies Rat significantly (2^, H3, while the heat- 
transfer channel remains virtually insensitive to the field. 
So, provided the intrinsic heating model is valid, one 
has to expect the heating-induced features of I-V, and 
the critical heat load in particular, to be field indepen- 
dent. The 'isothermal' I-Vs taken on the Bi2212 Dayem 
bridge at different fields (inset to Fig. 4 in is repro- 
duced in Fig.6A here) offer an independent test of the 
model. These data reveal a quasi-linear negative mag- 
netoresistance, whose value, SR sq (30T) / R sq (0) ~ —42%, 
correlates reasonably with direct measurements of the 
normal state MR in Bi2212 [US E3 Fig.6C shows 
the same data re-plotted as a sample resistance, R=V/I, 
normalised by its value at V— >0 versus the heat load, 
P=IV/A. With remarkable similarity to Fig.5B, these 
curves reveal a well defined threshold level P c , which 
separates the intrinsic and extrinsic contributions. More- 
over, as is clear from Fig.6C, different I-V merge onto the 
same curve, thus supporting the Joule heating origin of 
the deviations from I-V linearity. 

Moreover, one may note the qualitative similarity with 
the results of a similar analysis obtained for Bi2212 mesas 
[3l| . reproduced in Fig.6B. Here, the resistance of the 
sample is plotted together with the topmost electrode 
overheating above the base temperature, Tb- Fig.6B sug- 
gests the Ohmic response (the broken lines in Fig.6B) 
to occur in the absence of overheating. At higher heat 
loads, we see a progressive rise in overheating, which en- 
tails a corresponding reduction of Rn- The dotted line in 
Fig.6B represents the intrinsic 'heating-free' I(V), which 
has little in common with Josephson behaviour. The 
noticeable difference between the critical heat loads es- 
timated from the point of I-V departure from linearity 
for Bi2212 mesas and 3D-microbridges is probably due 
to differences in the heat transfer efficiency in these ex- 
periments. Besides, P c for the first case is most likely 
underestimated because of the unknown contact resis- 
tance. 
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FIG. 6: A The 'isothermal' IVC's of the Bi2212 Dayem bridge 
taken at at different fields from (3^1 : C The heat load depen- 
dence of the Dayem bridge' resistance, R=V/I, normalised by 
R N {B) for 10, 20, and 30Tesla (■,□,+ respectively). Insert 
B shows mutual correlation between the mean sample heating 
and deviations from I-V linearity from|3ll|. 



G. How does the heating-free I-V look? 

The results presented in Fig.5B,6B and especially the 
I-V scaling in Fig.6C provide a strong indication that the 
departures of experimental I-V from linearity originate 
in intrinsic heating. These results suggest the intrinsic 
('heating-free') out-of-plane response to be Ohmic. Such 
a conclusion drastically contradicts the naive picture of 
Josephson interplane tunnelling [32] but is supported by 
direct experiments [TtII^ performed sufficiently close to 
T c or the irreversibility field, Hi rr . 

In particular, comparison of the zero-field I-V mea- 
sured at the same base temperature on a single sample 
situated in a vacuum chamber and immersed in liquid ni- 
trogen reveals highly nonlinear Josephson-like branches 
in the first environment and practically linear ones in the 
latter case. I believe that the linear branches in this ex- 
periment are measured at an insignificant level of Joule 
heating of the structures because (i) the heat transfer co- 
efficient into liquid radically exceeds that into a vacuum 
and (ii) the characteristic heat load measured for these 
branches falls well below the threshold P c estimate. Re- 
markably similar results were obtained from the system- 
atic mixed-state I-V studies performed in the vicinity of 
Hi rr (T) at vastly different orientations of the magnetic 
field with respect to the (ab)-plane, as is illustrated in 
Fig.7. 

This similarity makes it possible to outline the distinc- 
tive features of the brush-like I-V measured at negligi- 
ble self-heating: (i) a voltage jump takes place from zero 
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FIG. 7: Typical I-V brush measured under conditions where 
the critical parameters of Bi2212 are suppressed by magnetic 
field ||(ab). Insets show corresponding differential resistances, 
7?#, and excess current, /*, of the branch number 



voltage at some characteristic current, Ij; (ii) all resistive 
branches have approximately the same excess current, /*, 
(given by the intersection of their slopes with the current 
axis); (iii) the differential resistance R# of a resistive 
branch is proportional to its number, and represents 
a fraction of the normal state resistance of the same junc- 
tion measured under conditions of complete suppression 
of its superconductivity. The remarkable similarity of 
these features with the characteristic manifestations of 
phase-slip centers of different dimensionality [l3L FLU ITo| 
in conventional superconductors suggests the phase-slip 
mechanism might be responsible for the 'intrinsic Joseph- 
son effect'. 

Evidently, the typical geometrical shape of a Bi2212 
sample is a thin flake rather than the quasi one- 
dimensional samples employed in classical PSC studies. 
However, those quasi- ID samples were made of almost 
isotropic material, in drastic contrast to Bi2212, which 
reveals extremely high and temperature-dependent elec- 
tric anisotropy, p c j p a b > 10 5 , [3(J. On taking account of 
this anisotropy, one immediately concludes that the typi- 
cal length-to- width ratio for an 'isotropic' Bi2212 sample 
falls into the 10-1000 range (i.e. is the same if not bigger 
than the typical ratios quoted in classical PSC studies 
EJ0Di so that our Bi2212 samples could be treated as 
quasi-lD. Nevertheless, Bi2212 seems to fall beyond the 
scope considered in [13 because the width of the sam- 
ple exceeds both the coherence length and the magnetic 
field penetration depth. However, the most recent exper- 
iment shed some light on the problem. The authors 
of 0| reported behavior very similar to PSCs (and to 
our findings) to occur in wide superconducting films (of 
the width d), in which both the coherence length and the 
magnetic field penetration depth are significantly smaller 



than d. In contrast with the ID case, it was found that 
oscillations of the order parameter may not necessarily be 
uniform along the two-dimensional analogue of a PSC. In 
particular, it was shown that these oscillations may occur 
in the form of propagation of waves carrying the order 
parameter singularities across the film. It seems plausible 
that likewise |l3l fLU Il5j | nonequilibrium superconductiv- 
ity of the phase-slip planes (PSP) might be responsible 
for the heating-free brush-like I-V reported here. 

In fact, the one-dimensional phase-slip center could 
naturally be considered as a subclass of weak links, char- 
acterised by direct conductivity (unlike that of well- 
known tunnel junctions). Such weak links have been ex- 
tensively studied, particularly for the practical realisation 
of Josephson electronics, cf. eg. review |T^. Contrary to 
tunnel junctions, weak links are characterised by a high 
current density, which generates a highly nonequilibrium 
state in the structure. A quantitative description of this 
state is rather complicated because it must address mu- 
tually correlated nonequilibria in three systems: the con- 
densate, the quasiparticles and the phonons of the struc- 
ture under consideration, together with its interaction 
with the environment. Despite the absence of any self- 
consistent theoretical description, the present-day level of 
phenomenological understanding |13l Il4l Il5l |26( appears 
sufficient to obtain a general idea about the underlying 
processes. Let us consider a model sample containing one 
(or many) PSC and account for the Joule heat generated 
in such a center. This heat can only escape through the 
surface of the sample by thermal conductivity along the 
superconducting material and/or directly into the sub- 
strate (or rather to the exchange gas). Heating is small 
when there is sufficiently good thermal contact and mod- 
erate Joule dissipation. In this case, the differential resis- 
tance of PSC is known to be determined by the normal 
state resistance, i?jv, |13|| which would be obtained on 
a single sample in conditions where its superconductiv- 
ity is destroyed somehow, for example by magnetic fields 
[lq j. In most cases, however, the thermal contact is far 
from ideal and this entails more or less noticeable heat- 
ing of the area near the PSC. If the dissipation in the 
phase-slip centre raises its temperature above the tran- 
sition temperature, a normal 'hot spot' 16] is created 
and an SNS junction is therefore formed. The Josephson 
effects are eventually quenched when the normal region 
grows. 

Surprisingly, this oversimplified scenario seems to be 
compatible with the experiment. Indeed, the differential 
resistance of the 'heating-free' branch correlates reason- 
ably with i?jv |10| . The increasingly strong deviations 
from linearity with the heat load correspond naturally 
to a case in which heating of the PSP becomes notice- 
able on the scale of Rn(T). As soon as the bias is suffi- 
ciently high, either the PSP or the whole sample becomes 
normal; however, its TV remain tunnel-like and reveal a 
clear 'gap' until the mean temperature of the stack ex- 
ceeds that of the minimum in R^(T), in agreement with 
[l8j . Finally, at exceptionally high levels of dissipation 
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and heating, the 'pseudo-gap' humps disappear and we 
obtain a slightly convex I-V, which is frequently identi- 
fied with Rjy. The mean temperature of the overheated 
stack, naturally, gives a very coarse estimate of the highly 
noncquilibrium temperature distribution across the stack 
and should therefore be treated with great caution 

III. CONCLUSIONS 

The intrinsic heating model proposed in is con- 
structed on the basis of Newton's Law of Cooling and 
the linearity of I-V experimentally obtained in the nor- 
mal state of the crystal. This model shows that the vast 
variety of results obtained with intrinsic tunnelling spec- 
troscopy originate from the Joule heating of the sample 
and are related to the temperature dependence of normal 
state resistance. In particular, the IT feature attributed 
to the normal state pseudogap in the electronic density 
of states is an indirect though robust signature of the 
true normal state R(T), while the 'superconducting' one 
is most likely caused by the 3-point configuration de- 
veloped in [3, E3, 13 • Normal state magnetoresistance 
is the most probable cause of the contrasting response 
to magnetic field of these two gap-like I-V peculiarities. 
Moreover, this scenario naturally explains why the 'su- 
perconducting' peak in 0, 0] collapses well below 
the superconducting critical temperature, T c , of material 



as well as why it survives a magnetic field which signifi- 
cantly exceeds the upper critical field of this substance; 
see 0, |22, 113, E3 f° r more detail. Thus, the conclu- 
sions gathered by temperature, magnetic field, intercala- 
tion and doping dependent IT are most likely irrelevant. 
Unfortunately, even the low-bias results reported by this 
group are of little help because of the 3-point configura- 
tion and unknown field-temperature-current dependence 
of the contact resistance. Indeed, analytical and experi- 
mental results suggest that, provided self-heating is neg- 
ligible, the individual branches of the brush-like I-V are 
linear. This finding does not support the conventional 
picture of 'intrinsic' Josephson tunnelling. Moreover, it 
seems plausible that, as in 0, Q, 0], stationary and 
nonstationary noncquilibrium superconductivity, phase- 
slip planes for example, might be responsible for the 
heating- free brush-like out-of-plane I-V in HTSC. 
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